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Abstract. We consider a diffuse interface model for an incompressible isothermal mixture of 
two viscous Newtonian fluids with different densities in a bounded domain in two or three space 
dimensions. The model is the nonlocal version of the one recently derived by Abels, Garcke 
and Griin and consists in a Navier-Stokes type system coupled with a convective nonlocal 
Gahn-Hilliard equation. The density of the mixture depends on an order parameter. For this 
nonlocal system we prove existence of global dissipative weak solutions for the case of singular 
double-well potentials and non degenerate mobilities. To this goal we devise an approach 
which is completely independent of the one employed by Abels, Depner and Garcke to establish 
existence of weak solutions for the local Abels et al. model. 
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1 Introduction 

In this paper we study the following nonlocal Gahn-Hilliard/Navier-Stokes type system 


{pu)t + div(pu (g) It) — 2div(z/((p)D'u) -|- Vtt -|- div('u ® J) = pVif -|- h, (1.1) 

div(u) = 0, (1.2) 

ift + u ■ Vif = div{m{ip)'Vp), (1.3) 

p = aLf - J * Lf + F'{ip), (1.4) 

J = -l3m{(p)Vp, f3 = {p2- Pi)/2, (1-5) 

P{^) = ^iP2 + Pi)+ ^iP2-pi)^, ( 1 - 6 ) 
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in Q := f2 X (0,T), where C h = 2,3, is a bounded smooth domain and T > 0 is an 
arbitrary hnal time. The associated boundary and initial conditions are 

= 0, = 0, on 912, 

on 

it(0) = ito, (^(0) = (^0, in 12, 

where 912 is the boundary of 12 and n is its outward unit normal. 

System fll.lD - fll.SI) couples a momentum balance equation fll.ll) for the velocity held u with 
a nonlocal convective Cahn-Hilliard equation fll.3p for the order parameter (p (difference of the 
volume fractions of the huids) and describes the how and phase separation of an isothermal 
mixture of two incompressible Newtonian viscous immiscible huids with diherent densities tak¬ 
ing into account long-range interactions between the molecules. Equation fll.2p accounts for 
the incompressibility of the mixture, pi,p 2 > 0 are the specihc constant mass densities of the 
unmixed huids, p = p{ip) given by 01.61) is the density of the mixture, tt is the pressure, h 
is the external volume force density and D denotes the symmetric gradient, which is dehned 
by Du := (Vn -|- V^rt)/2. Moreover, if a, 6 G M'^, we denote hj a ^ b the tensor dehned by 
(a (g) b)ij = ttibj, for i,j = l,---,d. 

The mobility m in 01.31) and the viscosity u in 01.11) are assumed to be </?—dependent and 
non degenerate, namely both are bounded from below (and above) by positive constants. The 
chemical potential p contains the spatial convolution J * p over 12, dehned by 

{J*p){x):= / J{x - y)p{y)dy, a; G 12, 

Jq 

of the order parameter p with a sufficiently smooth interaction kernel J satisfying J{z) = J{—z). 
Moreover, a is given by 

a{x) := / J{x- y)dy, 

Jq 

for x G 12. The double-well potential F is assumed to be singular and, in particular, a physically 
interesting case that will be included in our analysis is the following (see ini) 

^(f) = |((1 + s)>os(l+s) + (1 - - f)) - 0 < « < «„, (1.9) 

where 6, 6c are the (absolute) temperature and the critical temperature, respectively. 

System fll.lD - fll.bp represents the nonlocal version of the well known thermodynamically 
consistent dihuse interface model for two-phase how with diherent densities derived by Abels, 
Garcke and Griin in [8]. We recall that the local model deduced in [8] consists in the above 
system with the chemical potential p replaced with the local one 

11 = -A<f + (1.10) 


(1.7) 

( 1 . 8 ) 
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and completed with an additional homogeneous Neumann boundary condition for tp. 

We recall that the local chemical potential fll.lOp is the first variation of the local free energy 
functional (see HI) 

Actually, the local free energy considered in [8] contains also a positive coefficient a(<p) multi¬ 
plying the |V(pp under the integral. However, here we have set a{p>) = 1, since this coefficient 
does not introduce substantial complications into the analysis. 

A different form of the free energy can be associated to the ffuid mixture, more precisely 
the one proposed in [321EZ] and rigorously justified as a macroscopic limit of microscopic phase 
segregation models with particles conserving dynamics (see also [IS]). In this case the gradient 
term is replaced by a nonlocal spatial interaction integral, namely 

= \ [ [ J{x-y){ip{x)-(p{y)Ydxdy+ [ F{ip), 

^ Jn Jn Jn 

and the nonlocal chemical potential given by fll.dp is obtained by taking the first variation of E. 
The physical relevance of nonlocal interactions was already pointed out in the pioneering paper 
EZl (see also [Ml 4.2] and references therein) and studied (in the case of constant velocity) for 
different kind of evolution equations, mainly Cahn-Hilliard and phase held systems, see, e.g., 

[IDliinilSSllaSlETllSEllSSllSlIlElllglig 

Diffuse interface models for two-phase how of huids with identical densities are very well 
established and studied in literature. These models are based on the so-called model H (see 
[1311301, cf. also [ZBES] and references therein), in which the sharp interface separating the 
two huids is replaced by a dihuse one by introducing an order parameter (cf. [3]). They consist 
of the Navier-Stokes equations for the velocity held u nonlinearly coupled with a convective 
Cahn-Hilliard equation for an order parameter ip (cf., for instance, [HI SHI SIl 1131133113S113U]). 

As far as analytical results for the matched density case (i.e., pi = P 2 ) are concerned, 
the local Cahn-Hilliard/Navier-Stokes system has been tackled by several authors (see, e.g., 

[H 121131 [m HSl 13111321133113211321 ESI EHl El] and also [71 HH 133 137] for models with shear 

dependent viscosity), while the investigation of its nonlocal version (from the analytical view¬ 
point concerning well-posedness and long-term behavior) has started only more recently (cf., 
e.g., [13123123 EH 12S1E313D]). In particular, the following situations have been addressed: 
regular potential F associated with constant mobility in [13 123 123 ES]; singular potential as¬ 
sociated with constant mobility in [23; singular potential and degenerate mobility in [2H]; the 
case of nonconstant viscosity in [ 23 , which is particularly delicate as far as regularity results in 
two dimensions are concerned. In the two-dimensional case it was shown in [28] that for regular 
potentials and constant mobilities the problem fll.ip - fll.8l) with pi = p 2 admits a unique strong 
solution. Recently, uniqueness was proved also for weak solutions (see |2S])- Moreover, relying 
on the uniqueness results of [ 23 ] and | 23 ] a related optimal control problem was studied in [3D] 
for the case of constant mobility and regular potential. 
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Despite the considerable amount of contributions dealing with the matched density case, 
analytical results related to models for two-phase flow of fluids with unmatched densities are 
quite sporadic. In particular, as far as the local Abels-Garcke-Grun model is concerned, the 
hrst results on existence of weak solutions were obtained by Abels, Depner and Garcke in [5], 
for the system with singular potential and non degenerate mobility and in [B] , for the case of a 
regular potential and degenerate mobility. Regarding other diffuse interface models for fluids 
with different densities we recall the one considered by Boyer in [T3]. He proved existence of 
local in time strong solutions and existence of global weak solutions provided the densities of 
the fluids are sufficiently close. We also recall the quasi-incompressible model of Lowengrub and 
Truskinovsky [5S], where the velocity held is not divergence free, for which the hrst analytical 
results were obtained in [3111]. 

As far as nonlocal models for huids with unmatched densities are concerned, to the best of 
our knowledge no analytical results have been established so far and this paper aims to be a 
hrst contribution in this direction. More precisely, the goal of this paper is to prove existence 
of global dissipative weak solutions for the nonlocal Abels-Garcke-Grun model given by system 
(iriD-drsD. assuming, as in [5], that the potential is singular and the mobility is non degenerate. 
By weak solutions here we mean solutions with the minimum regularity requirement to allow 
a hnite energy and the validity of an energy dissipation inequality. 

Before explaining the strategy of the proof, let us briehy recall the approach used in [5] 
and discuss on the possibility to apply this approach to prove existence of weak solutions for 
the nonlocal system fll.ip - fll.Sp . In [5] existence of a weak solution is established by employing 
an implicit time discretization scheme. In particular, we point out that the Leray-Schauder 
hxed-point argument devised for the existence of a solution of the time-discrete problem (cf. 
[S] Lemma 4.3]) relies on the possibility of inverting the local relation between the chemical 
potential /i and if given by fll.lOp . This possibility is due to the fact that the relation between 
/i and f can be expressed by means of a maximal monotone operator since fi can be viewed as 
the subdifferential of the lower-semicontinuous convex (up to a quadratic perturbation) local 
functional Sioc- This approach allows in particular to keep f between the singular points —1 and 
1 in all the analysis. Indeed, the Abels-Garcke-Griin model is meaningful only when we have 
a bound on f between —1 and 1 which allows to keep the density p bounded from below and 
above by positive constants. This bound is ensured, in the case of the local Abels-Garcke-Griin 
model, by working with a singular potential as in |5] or with a degenerate mobility as in |6]. 

Now, the direct application of the approach devised in [S] seems hard in the present situa¬ 
tion. Indeed, the nonlocal chemical potential p can no longer be expressed as a sub differential 
of a lower semicontinuous convex functional and therefore the theory of maximal monotone 
operators is not directly applicable in the analysis. Moreover, the inversion of the nonlocal re¬ 
lation between p and f, under some reasonable conditions on the kernel J and on the potential 
F, seems to be a rather difficult task. 

A possibility to still exploit the approach of [S] in order to prove existence of weak solutions 
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for the nonlocal system could be to introduce a local perturbation term of the form —6A(p 
on the right hand side of fll.dp . Existence of a weak solution to the corresponding perturbed 
system could be proven, for every 5 > 0, by suitably adapting the argument of [5]. Then, 
existence of a weak solution to the original problem would be obtained by passing to the limit 
as 5 —>■ 0 (arguing as in Step III of the proof of the main result of the present paper). 

Although this approach would be possible, however, we propose here an alternative strat¬ 
egy which does not rely at all on the result of [S]. Our approach does not employ a time- 
discretization scheme and does not make use of Leray-Schauder hxed point arguments, but it 
is essentially based on the Faedo-Galerkin method and hence it is particularly suitable for a 
possible numerical implementation. 

Let us now describe the main lines of our approach. The starting idea consists in approx¬ 
imating the singular potential F by a suitable family of regular potentials dehned on the 
whole of M. This idea, that we already used in [27] for the same nonlocal system with matched 
desities, is quite classical (see, e.g., mm)- Nevertheless, it leads to some troubles when ap¬ 
plied to our problem. Indeed, if F is replaced by F^, we shall have to solve a problem in which 
the values of (the (p—component of the solution to the approximate problem with potential 
Fe) are no longer restricted to (—1,1) but belong to the whole of M. This implies that p{(pe) in 
this e—approximate problem is no longer a-priori bounded from below by a positive constant 
and consequently we are in trouble to get an estimate for the velocity held u^. 

To overcome this difficulty a possibility is to replace the linear density function p{(p) by a 
hxed smooth extension p{(p) from [—1,1] onto M satisfying 

0 < p* < p(s) < p*, < Rk, Vs e R, k = l,2, (1-11) 

p(s)=p(s). Vs e [-1,1], 

where p*,p*,Fi,F 2 are some given positive constants. However, we are now led to a further 
difficulty. Indeed, if we deduce (formally) an energy equation from system fll.lD - fll.Sp in which 
the linear function p is replaced by the nonlinear function p, by multiplying fll.ll) by u, 01.31) 
by p, integrating over H by parts and taking 01.4I) - 01.6I) . the incompressibility condition 01.21) 
and the boundary conditions into account, after some computations we obtain 

4 ( / + 2 f iy{p)\Du\‘^+ j m((p)|Vp|^ 

= 1 f ^'((p)m((p)(V(p-Vp)n2+ [ h -u, (1.12) 

^ Jn Jn 

where 

= T [ [ Jix-y){^ix)-p{y)Ydxdy+ [ F,((p), 

^ Jn Jn Jn 

and where we have denoted u^, simply by u, p, p, for the sake of simplicity. Therefore, a 
nonlinear p in system Ol.ip - Ql.Sp destroys the energy balance. A possibility to handle the nasty 
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nonlinear term on the right hand side of fll.l2j) is to recover the energy balance by inserting, in 
the approximate problem with the potential the term (l/2)p"(<yc)m((/?)(V<p • fi)u on the 
left hand side of the momentnm-balance equation fll.ip . This easy trick leads however to still 
another problem, namely, the problem to pass to the limit in this new “artihcial” nonlinear 
term. The idea at this point is to introduce some suitable regularizing terms in the system, 
depending on another positive parameter 6 which will be made go to zero in a second time. 
These regularizing terms, which allow to gain enough compactness to be able to pass to the 
limit, must be cleverly devised, since: (i) the energy balance should not be destroyed and (ii) 
when passing to the limit, hrstly as e ^ 0 and secondly as <5 —)■ 0, it should still be possible 
to prove that the limit (p satishes the bound |(p| < 1. This bound on ip will in particular 
allow the nasty artihcial term to vanish in the limit, permitting then to recover the original 
momentum-balance equation. More precisely, the regularizing terms that have been proven 
to be effective to our purpose are the term 5A^u in fll.ll) (this means, more exactly, that the 
term u, is introduced in the variational formulation of o with test function 

w E D{A^^'^), cf. Dehnition[21 here A is the Stokes operator with no-slip boundary condition), 
and still the term —6Aip in the expression of the chemical potential p. 

Summing up, our approach consists in proving existence of a weak solution to problem 
fll.ip - fll.6l) by approximating this problem with a two-parameter family of problems of the 
following form 


{pu)t -|- div(pw u) — 2div(i/((p)Zi)ri) -|- 6A^u -|- Vvr -|- div(ii (g) J) 


+ • V/x)^ = fiVip + h, 

( 1 . 13 ) 

dw{u) = 0, 

( 1 . 14 ) 

Lpt + u- Vcp = div(m(<yc)V/x), 

( 1 . 15 ) 

lj, = aip — J*ip + F^^tp) — 6Aip, 

( 1 . 16 ) 

J ■= -^((p)m((p)V/x, 

( 1 . 17 ) 

n = 0, ^ ^ = 0^ on T, 

on on 

( 1 . 18 ) 

'u(O) = Uq, V9(0) = p>os, 

( 1 . 19 ) 


where e and 6 are two hxed parameter. Notice that together with the regularizing term —6A(p 
introduced into the chemical potential, a homogeneous Neumann boundary condition for (p has 
to be introduced into the approximate problem and, moreover, the initial datum for ip has to 
be suitably approximated. The existence of a weak solution to the original problem will then 
be recovered by passing to the limit in two steps in fll.13p - fll.19p . i.e., by hrst passing to the 
limit as e —0 (with 6 hxed) and then as 5 —?■ 0. But, of course, before doing this we must prove 
that problem fll.13p - fll.19p (for e and 6 hxed) admits a weak solution. This will be achieved as 
hrst step by means of a Faedo-Galerkin procedure. 

The plan of the paper is as follows: in Section [2] we introduce some notation, recall some 
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classical results and preliminary lemmas. In Section |3] we formulate the assumptions, the 
dehnition of weak solution and we state the main result on existence of weak solutions. Section 
0] is entirely devoted to the proof of the main result. Since, as explained above, the proof is 
accomplished by a three level approximation of the original system. Section 0] has been split into 
three subsections for each step of the approximation argument: in Subsection 14.11 we develop 
the Faedo-Galerkin approximation scheme to prove existence of a solution to problem fll.l3j) - 
fll.lhp : in Subsection 14.21 we derive uniform in e estimates that allow to pass to the limit as 

e —)■ 0, in Subsection 14.31 we obtain uniform in S bounds, we shall pass to the limit as 5 —)■ 0 

and conclude the proof. 

2 Preliminaries 

Throughout the paper, we set H = V := and we denote by || • || and (•, •) the 

standard norm and the scalar product, respectively, in H as well as in L‘^{QY and 
The notations (•, ■)x and || • ||x will stand for the duality pairing between a Banach space X 
and its dual X\ and for the norm of X, respectively. 

We introduce the standard Hilbert spaces for the Navier-Stokes equations (see, e.g., [SH]) 

Gdi. := Vdi. := V := {n e C'o“(f2)" : div(n) = 0}, 

and recall that these spaces, for Lipschitz bounded domains, can be characterized in the fol¬ 

lowing way 

Gdiv ■= {u e L^{VtY : div(n) = 0, n • n\on = 0}, Vdiv := {n G Hq{QY ■ div(n) = 0}. 

The norm and scalar product in Gdiv will be denoted again by || • || and (•, •), respectively, and 
the space Vdiv is endowed with the scalar product 

{ui,U2)vdi„ ■= (V'Ui,Vn2) = 2{Dui,Du2), Vni,n2 G Vdiv 

We also introduce the Stokes operator A with no-slip boundary condition (see, e.g., [59]). Recall 
that A : D{A) C Gdiv Gdiv is dehned as A := —PA, with domain D{A) = H‘^{V,Y H Vdiv, 
where P ; P^Y^Y Gdiv is the Leray projector. Moreover, A~^ : Gdiv Gdiv is a selfadjoint 
compact operator in Gdiv Therefore, according to classical results, A possesses a sequence of 
eigenvalues with 0 < Ai < A 2 < • ■ ■ and Xj —)■ 00 , and a family {wj}j^jq C D{A) of 

associated eigenfunctions which is an orthonormal basis in Gdiv Moreover, by means of spectral 
theory the fractional operators A^ are dehned for every s G M with domains D{A^^‘^), which 
are Hilbert spaces endowed with their natural norm and scalar product. Recall that, since hi 
is assumed to be smooth, then we have D{A^P) ^ H^{V.Y, for all s > 0. 

We also recall Poincare’s inequality 

Ai ||w|p < ||Vn|p yueVdiv 
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We will also need to use the operator B ;= —A +1 with homogeneous Neumann boundary 
condition. It is well known that B : D{B) G H H is an unbounded linear operator in H 
with the domain 

D{B) = {</? e : dip/dn = 0 on 

and that B~^ : H ^ H is a. selfadjoint compact operator on H. By a classical spectral theorem 
there exist a sequence of eigenvalues pj with 0 < pi < /i 2 < • • • and /ij —)■ cxo, and a family of 
associated eigenfunctions Wj E D{B) such that Bwj = fij Wj for all j G N. The family {wj}j^^ 
forms an orthonormal basis in H and is also orthogonal in V and D{B). 

Furthermore, for every f E V' we denote by / the average of / over hi, i.e., / := |r2|“^(/, l)v 
(here |r2| stands for the Lebesgue measure of hi), and we introduce the spaces 

V,,= {vEV-.v = 0}, ■.= {feV':J= 0}. 

If m G C'(M) satishes m* < m{s) < m* for all s G M, with m^,m* > 0, then, for every 
measurable : hi —)■ M we can dehne the operator : V ^ V hj 

{B^u,v)v ■= / m{ip)Vu-Vv, 'iu,v eV. 

Jn 

For every measurable this operator maps V onto Vq and its restriction to Vq (still denoted by 
B^) maps Vo onto isomorphically. Let us denote by A/"^ : Vq Vq the inverse map dehned 
by 

= /, V/ G Vq and Af^B^u = u, Wu E Vq. 

As is well known, for every / G Vq and every measurable <p, Af^f is the unique solution with 
zero mean value of the Neumann problem 

—div(m(<p)VM) = /, in hi, 
g = 0, on dn. 

Furthermore, the following relations hold 

{B^u,Af^f)v = {f,u)v, yuEV, V/G^o', 

{f,Af^9)v = {g,Af^f)v = [ m{ip)ViAf^f)-ViAf^g), ^f,geVQ. 

Jn 

It is also easy to see that, for every measurable <p, we have 

VlfWv < WMJWv < ^II/IIk-, V/ e q, 

IIL II 

We end this section recalling threee lemmas that shall be helpful in the analysis. 

The hrst one is a simple lemma which will be useful for passing to the limit in the 
viscosity and mobility terms of the energy inequality. Its proof is left to the reader. 


( 2 . 1 ) 

( 2 . 2 ) 

( 2 . 3 ) 

variable 


Lemma 1. Let Q C N > 1, and let {/„} C L°°{Q) be a sequence such that ||/n||Lo°(Q) < C 
and fn^f strongly in L^(Q). Let {gn} C L‘^{Q) be another sequence such that gn ^ g weakly 
in L‘^{Q). Then fnQn fg weakly in L‘^{Q). 

The next lemma will be needed to prove the weak continuity of velocities with values in 
Gdiv If X is a Banach space, we denote by T]]X) the topological vector space of weakly 

continuous functions / : [0,T] —)■ X. 

Lemma 2. Let X, Y be two Banach spaces such that Y ^ X and X' Y' densely. Then 
L“(0, T; r) n ^([0, T]; X) ^ ^([0, T]-Y). 

The last lemma will be useful to deduce the energy inequality. 

Lemma 3. Let : [0, T) ^ M, 0 < T < oo, be a lower semicontinuous function and let 
P ; (0, T) ^ M be an integrable function. Assume that the inequality 

£’(0)ci;(0) + f £{T)u'{T)dT > f 'D(T)u{T)dT 

Jo Jo 

holds for all to G hT^’^(0,T) with oo{T) = 0 and cu > 0. Then, we have 

£{t) + J T>{T)dT < £{s), 

for almost all s G [0,T), including s = 0, and for all t G [s,T). 

For a proof of the last two lemmas see, e.g., [3]. 

Throughout the paper we shall denote by c, C,.. some nonnegative constants the value of 
which may possibly change even within the same line. Generally, the value of these constants 
depend on the parameters of the problem (e.g., F, J, z/, m, pi, G) and on the data Uq, ipQ, h. 
Further or particular dependencies will be specihed on occurrence. 

3 Main result 

In this section we state the main result on existence of weak solutions of system fll.ip - fll.8p . 
The assumptions on the kernel J, on the mobility m and on the viscosity u are the following 

(Al) J G IF^’^(M'^), J{x) = J{—x), a{x) := / J{x — y)dy>0, a.e. x G G. 

Jn 

(A2) m G and there exist m*,m* > 0 such that 

m* < m{s) < m*, Vs G M. 

(A3) u G and there exist z/*, i/* > 0 such that 

z/* < z/(s) < z/*. Vs G M. 
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As far as the singular potential F is concerned we shall work under the same assumptions as 
in [27]. More precisely, we assume that F can be written in the form 

A = + ^2, 

where Fi G 1,1), for some hxed integer p > 3, -F 2 G that the following 

conditions are satished 

(A4) There exist ci > 0 and eo > 0 such that 

> Cl, Vs G (—1, —1 + Co] U [1 — Co, 1). 

(A5) There exists cq > 0 such that, for each fc = 0,1, • • • ,p and each j = 0,1, • • • , (p — 2)/2, 

>0, VsG [l-eo,l), 

Ff'+')(s) > 0, Ff'+')(s) <0, Vs G (-1, -1 + eo]. 

(A6) There exists eo > 0 such that F^'^ is non-decreasing in [1 — eo, 1) and non-increasing in 

(—1, —1 + eo]. 

(A7) There exists Co > 0 such that 

F'\s) + a{x) > Co, Vs G (—1,1), a.a. x G hi. 

(A8) lim^^ii A[(s) = ±cx). 

Finally, the assumption on the external force h is 
(A9) /iGL2(0,T;f^jJ,forallT>0. 

Remark 1. Assumption J G Q^n be weakened. Indeed, the behavior of the kernel at 

inhnity is not essential. Alternative conditions are J G where := {z G : \z\ < 

(5} with 6 := diam(r2), or J G — hi), where hi — hi := { 2 ; G : 2 ; = x — p, x,p G hi} or 

also (see, e.g., [Tn] l 

sup/ (|J(x - p)|-P |VJ(x - p)|)dp < cx). 
xen Jq 

Remark 2. Assumptions (A4)-(A8) are satished in the case of the physically relevant loga¬ 
rithmic double-well potential 01.91) for every hxed integer p > 3. In particular, setting 

Fi{s) = ^((1 + s) log(l + s) + (1 - s) log(l - s)), F 2 (s) = 

then it is easy to check that (A7) is satished if and only if inffi a> 9^ — 9. However, note that 
other reasonable potentials satisfy the above assumptions (e.g., the ones which are unbounded 
at the endpoints). 
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Let us state now the notion of weak solution to Problem fll.ljl - fll.Sjl . 

Definition 1. Let Uq E Gdiv, V^o ^ with -F((^o) ^ and 0 < T < +oo he given. A 

couple [n, (^] is a weak solution to fll.ip - fll.Sp on [0,T] corresponding to [wojV^o] */ 


• u, (f and p. satisfy 

u E T]; Gdiv) n -L^(0, T; Vdiv), (3-1) 

if E L°°{0,T;H)nL‘^{0,T;V), (3.2) 

p = aip — J*ip + F'i'p) G L^(0, T;V), (3.3) 

{pu)t E LF^{0, T; D{Ay), ift e L2(0, T; V), (3.4) 

and 

if E L°°{Q), \ip{x,t)\ < 1 a.e. {x,t) E Q := kl x (0,T); (3.5) 


• for every ip eV, every w E D{A) and for almost any t E (0,T) we have 

{{pu)t, - {pu ® u, Dw) + {2u{ip)Du, Dw) - {u® J, Dw) 

= w) + {h, (3.6) 

{ift, fj)v + {m{ip)Vp, Vfj) = (u, ifVfj), (3.7) 

where J = —(dm{(p)'Vp E L^(0,T; H); 

• the initial conditions it(0) = uq, (^(0) = (fo hold. 

Remark 3. Notice that 03.211 and the second of 03.411 imply that tp E C'([0, T]; i7). Hence, 
thanks also to 03.ip . the initial conditions 'n(O) = no, 9?(0) = </?o make sense. 

We are now ready to state the main result. 

Theorem 1. Assume that (A1)-(A9) are satisfied for some fixed integer p > 3, and d = 2,3. 
Let no G Gdiv, P>o G L°°{Q) such that F{(po) E L^{Q) and |^ol < 1- Then, for every T > 0 
Problem Ol.ip - Ql.Sp admits a weak solution [n, yj] on [0,T] corresponding to Uq, (po such that 

ip E L°°{0,T-,LP{n)), (3.8) 


and satisfying the following energy ineguality 

S{u(t),ip(t)) + J (2\\yv{ip)DufdT +\\ym{ip)V,i\\'‘)dT < S{u{s),ip(s)) + j {h,u)v.,„dT, 

(3.9) 


for almost all s E [0,T), including s = 0, and for all t E [s,T], where 


£{u,^) := f ^p{ip)u^ + E{ip), 

E{p) ■■= J*ip)+ f F{ip) = ^ / f J{x 


y){ip{x) -ip{y)Ydxdy + 
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4 Proof of the main result 


The proof will be carried out in three steps. In the hrst step we shall consider the two parameters 
approximate problem ,5 given by fll.13p - fll.19p (with both e, <5 > 0 hxed) and shall implement 
a Faedo-Galerkin approximation scheme to prove existence of a global weak solution 
to satisfying an energy inequality. In the second step we shall consider only 5 > 0 hxed and 
deduce some uniform in e bounds for the weak solution (that we can now denote by [u^, ip^) to 
problem P^^^ which will allow to pass to the limit as e —?■ 0 in the weak formulation of P^^^ and 
to prove that the family of solutions converges to a solution (that now we can denote 

by [us, 9?^]) to Problem P^ given by 


{pu)t + div(p'u ig) u) — 2div(^i/((p)Du') + 5A?u + Vvr + div(n 0 J) 
div(ri) = 0, 


+ n ■ V<p = div(m(<p)V/i), 

p = aip — J*ip + F'{(p) — 6Aip, 

J := —(3m{ip)'Vp, 

- n ^ ^ 

^ ’ dn dn 

u(0) = uo, (p(0) = (pos, 


= 0, on dQ, 


pV(p + h, 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 


in which the potential F is singular. Finally, in the third step we shall deduce uniform in 6 
estimates for the family of solutions [ug, tps] to Problem P^ and pass to the limit as <5 —)■ 0 to 
prove that [u^^ips] converges to a solution to the original problem fll.lll - fll.Sp . 

In all the analysis we shall consider only the case d = 3. If d = 2 all the steps of the proof 
of Theorem [1] can be repeated (with strong convergences in stronger norms in comparison with 
the 3D case). However, the result of Theorem [T] does not improve substantially in 2D (see 
Remark ED . 


4.1 Step I. Faedo-Galerkin approximation scheme. 

For problem fll.1311 - 01.191) we shall consider the general situation of a regular potential F^, 
that in this subsection we denote simply by F, of arbitrary polynomial growth. Therefore, the 
assumptions we make for F are the following (cf. |17j i 

(RPl) F G Gq^(M) and there exists cq > 0 such that 

F"{s) + a{x) > Co, Vs G M, a.e. a; G D. 

(RP2) F G G^(R) and there exist Ci > 0, C2 > 0 and p > 3 such that 

F"{s) + a{x) > ci|s|^“^ — C2, Vs G M, a.e. a; G D. 
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(RP3) There exist C 3 > 0, C 4 > 0 and r G (1,2] such that 

< C3|F(s)| +C4, VseM. 

Remark 4. Since F is bounded from below, it is easy to see that (RP3) implies that F has 
polynomial growth of order r', where r' G [2, cx)) is the conjugate index to r. Namely, there 
exist C 5 > 0 and ce > 0 such that 

|^(s)| < csisr' + ce, VsgM. (4.8) 

Observe that assumption (RP3) is fulhlled by a potential of arbitrary polynomial growth. 

The assumptions on the kernel J and on the mobility m and on the viscosity v are the same 
as (Al), (A 2 ) and (A3), respectively. 

A weak solution to Problem is a pair [w, (pj satisfying fll.l3^ - 111.19p in the following 
sense 

Definition 2. Let Uq G Gdiv, L’o ^ L°°{il) with R(</?o) ^ T^(0) and 0 < T < +cx) be given. A 
pair [u^ip] is a weak solution to fll.13p - fll.19p on [0,T] corresponding to [ito, 9 ?o] if 


• u, ip and p satisfy 

ue^do.ri 1 Grdiy^ n L2(0,T;T)(R3/2)), (4.9) 

p G L°°(0, T; V) n L^(0, T; ^^(O)), (4.10) 

p = ap — J*p + F\p) — 5Ap G T^(0, T;V), (4-11) 

{pu)t G L"(0, T; D{A^/^y), pt e ^'(0, T; V), (4.12) 


for some k > 1; 

• for every ip eV , every w G D{A^i‘^) and for almost any t G (0,T) we have 

{{pu)t, - (pu (8) u, Dw) + {2u{p)Du, Dw) + 5{A^^‘^u, 

- {u®J,Dw) + {p)m{p){\/p ■ Vp)u,w) = -{pVp,w) + (4.13) 

{pt, 'iIj)v + {m{p)Vp, VV') = {u, pVfj), (4.14) 

where J = —J?{p)m{p)'Vp G L^{0,T; H); 

• the initial conditions ri(0) = Uq, (p(0) = Pq hold. 

We are now ready to state the following 
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Lemma 4. Let assumptions (Al)-(A3), (A9) and (RP1)-(RP3) he satisfied. Let Uq G Gdiv 
and ipQ eV such that F{(Pq) G L^{Ll). Then, for every T > 0 Problem admits a 

weak solution [u,(p] on [0,T] corresponding to such that 


^eL^{0,T-L^m, 


(4.15) 


satisfying the following energy ineguality 

+ E{^p{t)) + ^||V(/?(t)|p + 2^ \\^/iy{^p)DufdT + 6 \\A^/^ufdT 
+ < j^\p{To)ul +E{ipo)+ ^\\^Tof + j^{h,u)v,jijT, (4.16) 

for almost all t G (0,T), where we have set 

E{t) ■= J * p) + j E{p) = ^f j J{x-y){p{x)-p{y)fdxdy+ f E{p). 

Proof. Let us assume in addition that Pq G D{B). Existence of a weak solution as well as the 
energy inequality in the more general case oi po eV with E{po) E L^(f2) can be recovered by 
means of a density argument, in the same fashion as in m Proof of Theorem 1], by exploiting 
in particular the fact that, due to (RPl), F is a quadratic perturbation of a convex function. 
We introduce the family {'Wj}j>i of the eigenfunctions of the Stokes operator A as a Galerkin 
base in Vdiv and the family of the eigenfunctions of i? as a Galerkin base in V. We 

dehne the n—dimensional subspaces Wn ■= (wi, • • ■ , Wn) and := (-^i, • • • , Tfn) and consider 
the orthogonal projectors on these subspaces in Gdiv and JL, respectively, i.e., ;= Pvv„ and 

En '■= P>I'„- 

We then look for three functions of the form 


^n(t) = Tn(t) = ^ ^ (t)fij, Pn{t) = ^ (t)V'i, 

i=i i=i 


j=i 


that solve the following approximating problem 


{p{pn)Un)\wk) - {p{p 


''nj'^n ^ '^n 


„ Dwk) + 2 {u{pn)Dun, Dwk) + 


— / u 


L ‘ i^Jn ' 2 (P (pn)^(pn) (Vpn * V 


H“ 2 (pP') ^Pp) '^71') 


= -{(fn'^l^n.'^k) + P = 1, * * * ,n 

(Pn^V'/c) + (^(Pp)V/Xn, = (tinPn, 

Pp -Pp(^Pp J ^ (Pn p F (Ptt.) 


'p; '^k^ 

(4.17) 

I,--- ,n 

(4.18) 


(4.19) 
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(4.20) 

(4.21) 


Jn ■ P iSPn)'^iSPri)^Pn 
^n(O) V^On; '^n(O) 

where primes denote derivatives with respect to time, ipon = PnP>o, '*^n( 0 ) = PnUo, and hn G 
C([0, T]; snch that hn ^ h in L^(0, T; Moreover, we assnme that the fnnction p 

satisfying fll.lip is hxed snch that p G By writing 

{p{ipn)Un)' = p{^n)u'n + '^{pn)Unp'n 

and observing that fl4.18p can be written as 

= Pn{di-V{m{pn)^Pn) “ Un ' V<^n), 


it is not difficult to see that solving this approximating problem is equivalent to solving a 
system of ordinary differential equations in the 2n unknowns which can be reduced 

in normal form thanks to the fact that we have p{s) > p*, for all s G M, with p* > 0. Indeed, 
this condition ensures that, for every hxed n G N, the vectors '\/p((pn)mi, • • • , •\/p((p„)m„ are 
linearly independent and hence the Gram matrix {(p(<p„)mj, Wk^}j^k=i,..n, which appears in the 
hrst term on the left hand side of fl4.17p when it is written explicitly in terms of the unknowns 
aj^\ is not singular. 

By the Cauchy-Lipschitz theorem we know that there exists T* G (0, +oo] such that this 
system admits a unique maximal solution ;= := ,b^'^) on 

[0,T,:) with G C\[0,T*y,R^). 

We now multiply fl4.17p by fl4.18p by and sum over k = 1, ■ ■ ■ ,n taking fl4.19p and 
fl4.20p into account. In doing this we also observe that the following identity holds 

(^{p{pn)Uny,Un'^ = ^ J^^p{pn)ul + ^ J^p{pn)tul 
d f \ 1 /* 

= ^/ ^PMuI + - / ^{pn)Pn{diy{m{pn)^Pn) -Un-Vpn)ul. (4.22) 

Moreover, on account of the incompressibility condition div(ri„) = 0 and of the no-slip boundary 
condition = 0 on dQ, we have 


i^Pi}Pn)'^n ® '^ni D'^n^ / '^n ' {p^n ' V) ^pijpnj'^n) / ^ P(.V^n)') 

Jn Jn 

P I P^Pnj'^n ■ ('^n ' 

Jn 

and the last term on the right hand side of this identity can be written as 


u 


ut 


PiPnj'^n ■ „ p{}Pn)ij^n ' V) / o ^PiP'n)') 
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Therefore, we have 


_ f _ I f 

i^Pi.V^nj'^n ® / {u^ • V p(<y?^)) ~ / ("^^n ' ^ Pn)'^n' (4.23) 

./o ^ ^ ./o 


Furthermore we have 

- / Un-Qn- V)lt„ = Un- div(w„ (g) J) = / ■ [(divj^)^^^ + (J„ ■ V)Wn] 

Jo ./ o Jo 

1 


= / u„ div Jn + / Jn • V -w = / div J. 


n: 


/o 


1 

2 

1 

2 


div(p' (99„)m(v9„) 


'^\ipn)m{ipn)(y(Pn ' V/Xn) U 


2 _ F 

n 2 


((p„)div(m(v 9 n) V/rn)^t, 


(4.24) 


By means of fl4.22p - fl4.24p after some easy computations we then arrive at the following identity 


A 

dt 


]^p{ipn)ul + E{ipn) + ^||V(Pn||^) + 2 || 1 ^ + d\\A^^‘^Ur. 


+ ||\/m((p„)V/in|P = {hn,Un)vai.- 


(4.25) 


We can now integrate fl4.25p between 0 and t (after splitting the term on the right hand side on 
account of (A9)) and use (A2), (A3), fll.lll) . (RP2) and the fact that, since <po ^ D{B), then 
we have <pon Pq in L°°{Q), to deduce the following bounds 


1 ^n| L°°(0,r;Gdi„)nL2(0,T;yji„) < C, 

(4.26) 

'“n||L2(o,T;D(A3/2)) < C, 

(4.27) 

\\Pn\\L°°{0,T-,LP{Q)) A C, 

(4.28) 

A C, 

(4.29) 

1 Vpn| a C, 

(4.30) 

\\F{(pn)\\L<=°{0,T-,L^{n)) < C, 

(4.31) 


where henceforth in this proof we shall denote by C a positive constant such that 

C = (7(11^X011, ||<Po||v, ||-^(‘^o)||li(o), ||^||L 2 (o,r;y;'.j)- 

The constant C also depends on F, J, m*, p* and hi. Notice that these bounds hold at hrst 
instance with T = T*. However, since we have ||'U„(t)|| = |a("')(t)| and |l<Pn(^)|| = then 

the same bounds yield T* = +oo, i.e. T > 0 can be hxed arbitrary and the above bounds hold 
for every T > 0. 

Moreover, multiplying fl4.19p by —A<p„ in H and observing that —A(p„ belongs to the 
subspace we have 


(h'ri! A<p^) (Vp-T^, V<Pfi) (^QiCPn J * Pn T F (<Pri) ^A<p^, A(p^) 
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= (V<yCn, {a + F”{ipn))V(pn + a - VJ * Lfn) + (5|| A(/?„ 

> Co\\V^Pn\\^ - 2||VJ||Ll||V(^„||||(^n|| + 5\\/\^Pn\\‘^ 

> y II V(^n|P - c|l<^n|P + 5||A(^„|P, 

and therefore, combining this estimate with 

(V/i„, < ^11 V^Jnll^ + —II V/in||^ 

we get 

II V/i„|P > ||| - c||(^„|p + <5co||A(^„|p. (4.32) 

By employing 04.321) . and using a classical elliptic regularity result, from 04.30p and 04.28p we 
hence deduce the following estimates 


IIV^n||L2(o,T;y) A C, 

^^^^ll‘^n||L2(0,T;r/2(Q)) < C. 


(4.33) 

(4.34) 


As far as the control of the sequence of the averages {/U.„} is concerned, we have 




|(a<^„ - J * </?„ + F'{Lpn) - 5ALpn, 1)1 = |(F'((^„), 1)1 


JQ 

< c||F(<^„)||ii(Q) +c<C, 

due to the bound 04.311) . By means of Poincare-Wirtinger inequality, this control and 04.301) 
imply the bound 


||hn||L2(0,T;y) A C. (4.35) 

Therefore, since 5 > 0 is hxed, from the estimates obtained above we deduce that there exist 
u, if and fi such that, up to a subsequence we have 



U, 

weakly* in L°°(0, T; Gd*^), 

(4.36) 



weakly in L^(0, T; Z1(A^'^^)), 

(4.37) 



weakly* in L“(0, T;V H L^(fl)), 

(4.38) 



weakly in L^(0, T; i7^(fl)). 

(4.39) 

l^n 


weakly in L^(0,T;1/). 

(4.40) 


As next step we need to derive some estimates for the two sequences of time derivatives 
{dtPn{p{^n) Unj } and Let us begin with the hrst sequence. Take w G D{A^P) and write 

w = Wj + Wjj, where Wj G and Wjj G (recall that Wj and Wjj are orthogonal in all 
Hilbert spaces D{A^P) for all 0 < s < 3). From 04.17p we can write 

(^dtPn{p{p>n)Un),W^^^^^^^^ = (p((p„)nn)', 
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= [p{Lpn)Un ® Un, DWi) - 2{l>{pn)DUn, DWj) - 6{A^''‘^Un, 


Un ■ {Jn ■ V)lt>/ - -(p"(v9„)m(v9n)(Vv3n ' VPn)Un, Wj) 
) ^ 


2 ^ V^n) ^V^n) '^ni '^1^ 

- ^(^^{^n){div{m{ipn)Vpn) - Pn{dw{m{ipn)VPn)))Un, Wj^ 

- ((^nVyUn, Wj) + (/l„, Wi)v^,^. 

We now estimate individually the terms on the right hand side of fl4.4ip . We have 


\[p{iPn)Un®Un,DWi) \ < c||||^||lU/1|^3(0)3 < c| 


u„\\ \\w 


ll_D(A3/2)) 


2\{u{^pn)DUr„Dwj) \ < c||V'u„||||m/||j^i(f^)3 < c||V'u„||||m||^(^i/ 2 ), 
S\{A^^‘^Un,A^^‘^Wj)\ < 1110/11^(^3/2) < (5||A^/^n„||||m||//(^3/2), 


Un ■ {Jn ■ V)lO/ 


Ur. 


{^{Pn) m{p>n)Vpn ■ '^)Wl 


< c||n„||||V/in||||lt>/||//3(f2)3 < c\\Un\\\\VPn\\\\w\\D{AV^), 

\{'^'{Pn)'m{ipn){'^^n ' Vpn)Un,Wj) < c|| V<^n ||L10/3(q)3 || V/i„ || || ||L6(q)3 || 10 1|//(A), 


(^^{^Pn){Pn{Un ' Vf^n) “ Un ' ViPn)Un,Wi^ < c\\Un ' Vf^n || || || H/1| H2 (q)3 


— c|| II /^6(/2)3 II VII L3(f/)3 II "lAfi II II ro II /)(^), 

^ ^p'(<^n)(div(m(<^„)V/i„) - P„(div(m(<^„)V/i„)))it„, 10/ 


(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 


< \\dw{m{pn)'^l^n)\\y,\\'^{Vn)Un ' U)i\ 

< c||V/i„|| (||lt„||i6(Q)3||lO/||//3(Q)3 + II Vn„|| ||tO/||//2(f/)3 + ||V<yC„||/^10/3(Q)3||n„||/^6(Q)3||?O/||//2(Q)3) 

< c||V/r„||(||VWn|| + ||V(//„||il0/3(f/)3||Wn||L6(O)3)||w7||/)(^3/2), (4.48) 

|(v?„V/i„,'fO/)| < |l9?„||||V/i„||||m/||//2(f^)3 < ||9?„||||V/i„||||'fo||//(A), (4.49) 

where fll.lip . (A2) and (A3) have been used. We now need the following interpolation embed¬ 
dings 

L°°(0, T; L\n)) n L\0, T; H^Q)) ^ L"(0, T; P®/"(fi)) L"(0, T; 

for 4 < s < cx), and 

L°°(0, T; H\n)) n L\0, T; H^Q)) ^ L"(0, T; H^+'^/^{n)) ^ L"(0, T; 
for 4 < s < oo. In particular we have 

L°°(0, T; L2(f2)) n L2(0, T; H%n)) ^ L%Q), L°°(0, T; H\n)) n L\0, T; P^(fi)) L^°{Q). 


\v 
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These interpolaton embeddings and fl4.26jl - fl4.30jl . fl4.34p entail the following bounds 


(4.50) 

(4.51) 


||'*^n||L6((3)3 < IIVWn||Ll8/5(Q)3x3 ^ 

|l‘^n||l,lO(Q) < Cs, l|Vv3n||LlO/3(Q)3 < C 5 , 

where the second bound in fl4.50ll follows from the bound of Vun in L®(0, T; m- 

L^(0, T; (for s > 12/5; take s = 18/5), and the second bound in fl4.5ip follows 

from the bound of V(pn in T^(0, T; if^/'^(r2)^) L^(0, T; (for s > 4/3; take 

s = 10/3). Therefore, by means of estimates fl4.42p - fl4.49p and on account of the hrst bound 
fl4.50p and of the second bound fl4.5ip . from fl4.4ip we then deduce the following estimate (not 
uniform in 5) 

||^t.Pn(p(^n)'*^n) ||i30/29('Q 7’.£)(^3/2y) < c*,. (4.52) 

Now, we have 

V [p{ipn)Un) = p{ipn)VUn + /S'((p„) V<Pn ' 

and hence from fl4.50p . fl4.5ip we see that V(J){(pn)Un) is bounded in L^^/’^(Q)^ which implies 
that 


||p(<yCn)u„||^i5/7(0,-r.,yi,i5/7(Q)3) < Cs. (4.53) 

To get strong convergence for the sequence of </?„, let us hrst observe that, taking '0 G 14 and 
writing + xpn, where -0/ ^ lE'n and 'll:!! G T/; (recall that -0/ and -0// are orthogonal in 

all Hilbert spaces 11(5^/^), for 0 < r < 1), from fl4.18p we have 

{^n^'^)v = (<^ni'0/) = -(V/i„, V?//) + {pnUn,ViJl), 


and, since p > 3 

\{(PnUn,\/i>i)\ < ||<^n||LP(0)||^in||L6(0)3|lVV’/|| < c|| (^„ || i,p(o) || || ||?/1| y, 

whence we deduce that 


|l‘^nlU2(0,T;V') < C. (4.54) 

Since pn is bounded in, e.g., L^(0, T; iJ^(r2)), by Aubin-Lions lemma we then deduce that 

—)■ (p, strongly in L^(0, T; iJ^“'^(r2)), 7>0 (4.55) 

and in particular we have 


Pn ‘P: pointwise a.e. in Q. 
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Observe also that by Lebesgue’s theorem we have 

p{^n) p{^), strongly in L''(Q), Vg G [2, cx)), (4.56) 

and the same strong convergence holds also for to l/p^ip). 

We now derive strong convergence for the seqnence of To this aim notice first that from 
fl4.53p we have 

ll-Pn(p(v^n)'*^n) ||2^2(o,r;irl(0)3) —^<5- 

From fl4.52p and fl4.57p . again by means of Anbin-Lions lemma we therefore dednce that 

Pn{p{^n)Un) ^ p{^)u, strongly in L^{Qf, 

for some p{(p)u G L‘^{Q)^. Bnt, from fl4.56p and the weak convergence available for Un we 
have p{pn)un p{p)u, weakly in L‘^{QY (also in L®“'^(Q)^, for 7 > 0 ), which implies that 
Pn{p{pn)Un) ^ p{.v)'i^-i Weakly in L‘^{Q)^. Therefore, we dednce that p{p)u = p{p)u and so 

Pn{p{^n)Un) strongly in L'^{Qf. (4.58) 

In particular, from fl4.52p . fl4.58p and from fl4.54p . fl4.38p there follows that up to a subsequence 
we have 


dtPn{p{^n)un) ^ {p{p>)u) Weakly in (0, T; T*(A^/^)'), (4.59) 

pti weakly in L^(0,T; V'). (4.60) 

With fl4.58p and strong convergence for the sequence of (pn at disposal we can now establish 
strong convergence for the sequence of by using a classical argument (cf. pQ, [511 Section 
2.1]), which we report for the reader’s convenience. Indeed, we have 

[ [ P{^n)ul = [ [ Pn{p{p>n)Un) ■ Un [ [ P{pW , (4.61) 

Jo Jn Jo Jn Jo Jn 

which means that the L^(( 5 )^-norm of ^J'p{ipn)un converges to the L^(Q)^-norm of ^J'p{ip)u. 
Since, thanks to strong convergence for the sequence of we have also ^Jp{pn)Un \/p{Pn)u, 
weakly in L‘^{Q)^, then \/p{ipn)un —t \/p{p>n)u, strongly in L‘^{Q)^. Finally, notice that due 
to fl4.56p we have 'p{ppn)~^P —f 'p{pp')~^P^ strongly in L^{Q) (for all q G [2, oo)) and therefore 
we get Un ^ u strongly in (7 > 0 ), which implies that, up to a subsequence we have 

Un ^ u pointwise a.e. in Q. Since the sequence of Un is bounded in L^{QY, then 

Un —)■ u, strongly in (4.62) 

Furthermore, from fl4.55p and the second fl4.5ip we obtain 

Wpn —t V</?, strongly in (4.63) 
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Now, the strong convergences fl4.62jl . fl4.63p . fl4.56p . together with the weak convergences 04.36^ - 
fl4.40p . fl4.59p . fl4.60p and with the strong convergences and p'((pn) —t in 

L‘^{Q), for all g < oo allow to pass to the limit in the approximate problem (14.1711 - 04.2ip and 
to recover the weak formulation 04.13p . 04.1411 . In particular, observe that we have 

Un ■ V(p„ —u ■ V(p, strongly in L‘^{Q), 

(also strongly in for all 7 > 0, due to 04.62p . 04.6311 1 and therefore 

Pn{un ■ V(pn) ^ u ■ V(p, strougly in L^(Q). 

Hence, when we pass to the limit the contribution of the seventh term on the left hand side of 
04.17P converges to zero. Moreover, we have 

div(m((p„) V/in) ^ div(m(y9) V/i), weakly in L^(0,T; V'), 

and therefore also 

Pn{div{m{ipn)VPn)) ^ div(m(<y9)V/i), weakly in L^(0,T;I/'). 

In addition, for all w E it is easy to see that we have 

{Lpn)Un'W -> {^p)uw, strougly in L^(0, T; V) 

(also strongly in T; due to 04.62p . 04.63p h and therefore also the con¬ 

tribution of the last term on the left hand side of fl4.17p converges to zero when passing to the 

limit. 

We now claim that u E Cu,([0, T]; To prove this claim, hrst observe that from (14.2611 

and (I4.52P we deduce the boundedness of 

Pn{p{‘^n)Un) iu L°°(0, T; Gd^) and 

Pn{p{^n)un) in D{A)') G([0,T];D(H)'). 

Therefore, thanks to Lemma [2] and on account of (14.5811 we infer that 

z ■.= p{p>)u E C^{[d,T]-L\nf). 

Moreover, since E G([0,T];iL), then, by (II.lip we have p(<p) E G([0,T];iL). Now, let u be 
the representative, in the equivalence class of n, given by n = p{(p)~^z. Due to the boundedness 
of p(<p) and noting in particular that we have also p(<p)“^ E G([0, T];iL), it is now easy to get 
the desired claim. 

Let us now prove that the initial conditions (11.1911 are satished. The argument to see that 
<p(0) = (fo holds, by integrating (I4.18p in time between 0 and t and using the strong convergence 
—)■ in G([0,T];iL) (which follows from (14.2911 and (I4.54p l is quite standard and we omit 


21 




















































the details. We just give some details on the argument to prove that n(0) = Uq. Let us then 
take w G and set := 'Yl!k=i^kWk-, for N < n, where = {w,Wk)- We multiply 

(14.17p by ak, sum on k from 1 to < n, and then integrate the resulting identity between 0 
and t to get 


p{ipn{t))Un{t),W^] - (p{ip0n)U0n,W^) - / {p{(Pn)Un ® Dw^)dT 


+ 2 / {y{!pn)DUn^Dw^)dT + 6 / {A^^'^Un, A^^^W^)dT 



Ur. 


n-{Jn-V)w^ + - / (p"(<y5„)m(<y5„)(V9?n • V/i„)n„,m^)cir 

0 Jn ^ Jo 

^ / (^{(Pn){Pn{Un-'^^n) -Un-V(Pn)Un,W^^dT 


{div{m{ipn)Vpn) “ Pn(div(m((p„) V/Xn))) it„, ] dr 


foi; 

= -/ (<yCnV/i„,m^)dr + / {hn,w^)dT 

Jo Jo 


(4.64) 


Let us now multiply (14.641) by y G C^(0,T) and integrate again in time between 0 and T. 
We then pass to the limit in the resulting identity hrst as n —?■ oo, by using (I4.58p and the 
weak/strong convergences obtained above, and then as iV —)■ oo. We then perform a similar 
computation on the weak formulation (I4.13p (with the same test function w) by hrst integrating 
it in time between 0 and t, by multiplying then the resulting identity by y and integrating again 
in time between 0 and T . Comparing the two identities obtained in this way we are led to 


p{ipo)uo,w 


X{t)dt = (p{lp{0))u{0),w) [ xit)dt, 


which yields p(</?( 0 ))n( 0 ) = ]){ipq)uq and therefore w( 0 ) = uq, since </?( 0 ) = ipQ. 

We now want to show the energy inequality (I4.16p . To this aim we integrate (14.251) in time 
between 0 and t and get 


[ ^Pi^n{t))ul{t) + E{(pn{t)) + ^\\ViPn{t)f + j (^2|| y/ p{ipn)DUnf + 5\\A^PUr,\\^ 

+ \\^Jm{ipn)y Pn\f)dT = J ^p{ipon)ulr, + E{p)0n) + {hn, Un)vdi^dT, (4.65) 

for all t G [0,T]. Then, we pass to the limit in (I4.65P taking the weak/strong convergences 
above into account and using the weak lower semicontinuity of norms. In particular, as far as 
the term containing the variable viscosity is concerned, we observe that Du^ Du weakly 
in while \/y{i~Pn) strongly in for all 2 < g < oo. Hence 

^ v[iPn)Dun \/v{i^Du Weakly in (7 > 0) and therefore also weakly in 
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By also employing Lemma[T]to pass to the liminf in the term containing the variable 
mobility we deduce 


2\\\/i'{(p)Du\\'^dT < liminf / 2|| \/i'{(pn)DUn\\‘^dT, 

n^oo Jq 

rt 

\\^/m{ip)VnW^dr < liminf / ||y/m((^„)V/in||^dr. 


Concerning the hrst term on the left hand side of fl4.65p . in view of fl4.62p and fl4.56p which 
imply that 


pi^n)ul p{(p)u‘^, strongly in 


(4.66) 


we have that this term converges for almost all t G (0, T) to the hrst term on the left hand side 
of fl4.16p . The fact that the term in the norm of V(pn converges, for almost all t G (0,T), 
to the third term on the left hand side of 04.161) is a consequence of 04.63p . The passage to 
the limit in the other terms in 04.65p is straightforward and this concludes the proof of 04.16p . 
which holds for almost all t G (0,T). 

Finally, we deduce an auxiliary energy inequality which shall turn out to be useful to deduce 
the energy inequality 03.9p . Let us multiply 04.25P by a; G hF^’^(0, T), with u)(T) = 0 and a; > 0. 
We obtain 


2pi^0n)ul^ + E{LpQn) + ^ )^( 0 ) + 


+ E((pJ + ^-\\Vpnf](^'{T)dT 


2|| \/i'{^n)DUn\\^ + + II v/m((^„)V/in|Pju;(r)dr - j (h,„, Un)vii^^{T)dT. 

(4.67) 

In order to pass to the limit in the term containing the functional E{(pn) on the left hand side 
of 04.67P we need to prove that 


/ u'{T)dT / E{ipn) / u'{T)dT / E{ip). (4.68) 

Jo Jn Jo Jn 

To this aim observe that, thanks to the following compact and continuous embeddings 

T; H\n)) n T; V') ^^(0, T; 77^(12)) 7:^(0, T; C(n)), 3/2 < s < 2, 

and to the bounds 04.341) and 04.54p we have 

^nir) (p{t), in C(f2), 

for almost all r G (0,T). On the other hand, due to the energy identity 04.65P the sequence of 
integrals E((pn(T)) is uniformly bounded with respect to n G N and for a.a. r G (0,T) (cf. 
04.3ip ). Hence, 04.68P follows immediately by applying Lebesgue’s theorem. 
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Passing to the limit in fl4.67jl and employing weak/strong convergences for Un, ^n, hn (recall, 
in particular, fl4.63p and fl4.66p h and weak lower semi continuity of norms in the same fashion 
as done for the proof of 04.1611 we get 

( f + E{^o) + ^||V(poir)^^(0) + f + E{^) + ^||V(pf)u;'(r)dr 

«/ w 0 V 

T T 

> {2\\^/v{ip)Du\\^ + 5\\A^/‘^u\\‘^ + ||^m(<p)V/i||^ja;(r)cir - (/i, n)v^,„a;(r)dr, (4.69) 

for all uj e iy^’^(0,T), with oj{T) = 0 and cn > 0. 

□ 


4.2 Step II. Limit as e —0. 

Next, we consider problem 04.ip - 04.7p where now the potential F is singular and 5 > 0 is still 
hxed (in this subsection we shall denote the initial datum for (p simply by (po, instead of (pots)- 
We aim to prove that this problem admits a weak solution by approximating it with a sequence 
of problems of the form 01.1311 - 01.1911 with regular potentials F^. More precisely, we prove 
the following 

Lemma 5. Let assumptions (Al)-(A9) he satisfied for some fixed integer p > 3. Let Uq G Gdiv, 
ipo eV nL°°(r2) such that Fi^ipo) G L^(r2) and |^ol < 1- Then, for every T > 0 Problem 04.ip - 
04.7p admits a weak solution [u,(p] on [0,T] corresponding to [uq,p)q] such that 

t.eC.(|0,T] 

; Gdiv ) nL2(0,T;D(7l3/2)), (4.70) 

(p G L°°(o, T-vn L^in)) n ^^(o, t- i7^(fi)), (4.71) 

peL\t),T-V), (4.72) 

satisfying the bound 

|<p(a:,t)| < 1, for a.e. {x,t) ^ Q, (4.73) 

and satisfying the energy ineguality 04.16p (with p in place ofp) for almost all t G (0,T). 

Proof. We consider Problem P^^^ consisting in 01.13p - 01.19p where the regular potential F is 
taken equal to F^ given by (see [27] and [29] i 

Fe = Fi^ + F2e, 

where Fi^ and F 2 e are dehned by 

— el s > 1 — e 

d'h), kl < 1 - € (4.74) 

fP(-1 + 6). s < -1 + t. 
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(4.75) 


( F^il-e), s>l-e 
= < F”{s), |s|<l-e 

[ + s <-1 + e, 

and Fi,(0) = Fi(0), F(,(0) = F((0 ),.., = ^(^”'^(0), and ^ 2 ,( 0 ) = ^ 2 ( 0 ), ^^(0) = 

^ 2 ( 0 ). Recalling [271 Lemmal, Lemma 2] and [221 Proof of Theorem 2], there exist two constants 
Cp and Dp, depending on p but independent of e, and there exists cq > 0 such that 

Fe{s) > Cp\s\P - Dp, Vs e M, Ve e (0 , cq], (4.76) 


and 


F"(s) + a{x) > Co, Vs G M, a.e. x G 12, Ve G (0,eo]. (4.77) 


Now, thanks to Lemma [4] we know that Problem P^^^ admits a weak solution [u^,Pf\ with the 
regularity properties fl4.9p - fl4.lip satisfying the energy inequality fl4.16p 


j ]^p{p^)ul ++ ^-\\Vipe\\‘^+ 2 j \\^Jv{^p^)Du^\\‘^dT + 5 f 

V «/ 0 *^0 

+ \\^/m{ip,)Vdr < j^^p{po)ul + E,{(po) + ^||Vy^oll^ + ^ {h,u,)v^ijT, 

for almost all t G (0,T), where 


(4.78) 


Fe{^) ■=^\\Va^\\‘^*(p) + [ E,{(p) 

tJ 

J{x-y){ip{x)-ip{y)Ydxdy+ f E^{ip). 



O JQ 


Notice that in the weak formulation of fll.ldp we don’t have additional auxiliary terms, like the 
last two terms on the left hand side of fl4.17p . 

Since 5 > 0 is hxed, starting from fl4.78p . splitting the last term on the right hand side by 
taking (A9) into account, and arguing as in the Faedo-Galerkin scheme of the proof of Lemma 01 
by employing assumptions (Al)“(A3), fll.lip . fl4.76p . fl4.77p and (A6) which implies that there 
exists Co > 0 such that 


FUs) < Fi(s), Vs G (-1,1), Ve G (0,eo], (4.79) 


we still deduce estimates fl4.26p - fl4.30p and fl4.33p . fl4.34p . fl4.54p for u^, and /i^, namely 

\\Ue\\L^{0,T-,Gaiv)nL^{0,T-,Vaiv) ^ F, (4.80) 

'“ell 1,2(0,T;D(A3/2)) ^ C, (4-81) 

IIVe||l,°°( 0 ,T;I,P(O))nL 2 ( 0 ,T;V) < C, (4.82) 

^^^^IIVe||i:,°°(0,T;V)nL2(0,T;r/2(Q)) < C, (4.83) 
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II V/ie|lL2(o,r;H) < C, 

\We\\L'^(0,T-,V') < C, 


(4.84) 

(4.85) 


which are now uniform in e. Here, all constants C have the same kind of dependencies on the 
data wq; '^Oi h and on the parameters of the problem as in Step I. Moreover, the estimate for 
the time derivatives ^ which corresponds to fl4.52p and obtained by comparison in the 

weak formulation of fll.ldp (cf. fl4.42p - fl4.49p i now becomes 

II j;,30/29(o,T;D(A3/2)') —^5- (4.86) 

We also need an estimate for the sequence of and in particular we need to control the 
sequence of averages /I^. To this aim we notice that equation fll.lbp can be written in the form 

(p( + li, • V(p, = (4.87) 

Test fl4.87p by . On account of fl2.1l) and fl2.2p we obtain 

(F,'(^,) - + {F'M - FlM.KXn. ■ VVe))y 

= (4.88) 

and since the elements in all dualities of this identity belong to 77 (the fact that F'^{ip^) G H 
follows from comparison in the expression for /i^ = aip^ — J + ~ then in fl4.88p 

we can replace all dualities in V with scalar products in H. Now we have 

- T;'(<Pe), = {Ki.Ve) “ a(p, - J * + F'^ifi^) - F/((p,) - (5A(p,) 

F-\\FlW,)-liM\X-hav,-J*V,f + S f F.'V,)|V^y 

«/ 

> - Cj\M^ + co6\\Vip4\^ - aoc<5||V<^,|ll (4.89) 

Therefore, by combining fl4.88l) with fl4.89l) we deduce 

WKiFe) - F’{(p,)\\ < c(||A/'^,<^(|l + ||A/'^,(ne • V<^,)|| + V6\\(p,\\v) 

< c(||<yc(|ly' + \\u^ ■ VipeWv + V^||<Pe||v), (4.90) 

where fl2.3p has been taken into account as well. Notice also that we have 
\{u, • V(^e,'0)l = l(we<Po V^/^)| < c||Vw,|1||(/?,||lp(o)||'0||v, 
for all 'll) & V, which yields, on account of fl4.80p and fl4.82p 

ll'l^e • V9?e|lL2(o_r;V') < C. (4.91) 

Hence, due to fl4.85lh fl4.9ip and fl4.82l) from fl4.90l) we obtain the uniform in 5 and e bound 

WKi-Ft) - -^e(<Pe)||L2(0,T;H) < C. (4.92) 
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With fl4.92p available, by employing the condition |<yCo| < 1 and the bound 


< \F[is)l ys e (-1,1), V6 e (0,eo], (4.93) 

for some eo > 0, which is ensured by (AS), (A6) and (A8) (see [271 Proof of Theorem 1] for 
details), we can now apply an argument devised by Kenmochi et ah |15] (see also [18] and m 
Proof of Theorem 1]) and deduce the following control 

\\K{ye)\\L2(0,T-,Li{n)) << 7 (^ 0 )- (4.94) 


Since 


He= {a(pe - J * - 6A(f,) = / F^'(<^e), 


then we have ||/rg||L2(o,r) < C and therefore, by Poincare-Wirtinger inequality we get 


||he||L2(0,T;V) < C. (4.95) 

Like in the Faedo-Galerkin scheme (5 > 0 here is still hxed), from estimates fl4.80p - fl4.84p and 
fl4.95p we deduce that there exist it, if and fi such that, up to a subsequence, we have 

^ u, weakly* in L°°(0, T; Grfj^), (4.96) 

ite ^ n, weakly in L^{0,T-, (4.97) 

f, f, weakly* in L°°(0, T;V f] L^(fl)), (4.98) 

weakly in L^(0, T; 77^(12)), (4.99) 

yUe ^/i, weakly in L^(0, T; 1/). (4.100) 

Furthermore, by using the bound fl4.52p and the bound (cf. fl4.53p l 

l|p(<7c)'“e||Ll5/7(0,r;Wl'l5/T(Q)3) < Cs, (4.101) 


as well as the bound 04.851) and arguing as in the Faedo-Galerkin scheme of Step I we can again 
deduce the strong convergences 


ife ^ if, strongly in L^(0, T; 77^ ^(^))) 7 > 0) (4.102) 

ite —)■ u, strongly in (4.103) 

Vife Vf, Strongly in L^°/3-^(Q)^ (4.104) 

as well as 

p"(</?,) p"(</?), rn{f,) m{f), strongly in L'^{Q), Vg G [2, cx)). (4.105) 


Notice that now, in order to deduce strong convergence for (and therefore fl4.103p and 
04.1041) 1 we do not need to employ the trick used in Step I (see 04.611) 1. In particular, fl4.86p . 
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fl4.101jl and Aubin-Lions lemma entail strong convergence for in (notice that 

^i,i5/7(q) 3 and that ^ D{A^/^y, since 

In order to pass to the limit in the variational formulation of Problem and hence to 
prove that [it, p] is a weak solution to Problem fl4.ip - fl4.7p we need to show that the limit ip 
satisfies the condition \ip\ < 1 a.e. in Q = hi x (0,T). This can be done exactly as in [271 
Proof of Theorem 1] by adapting an argument devised in [20] (see also [23]). We recall that 
this argument is based only on the use of: (i) estimate ||F^'((pe)|| 2 ,i(Q) < Cippo) (cf. fl4.94p h (ii) 
the pointwise convergence ip^^ ^ ip a.e. in Q (cf. fl4.102p l and (iii) the fact that F\s) ±cxo, 
as s ^ ±1 (cf. (A8)). Hence 04.731) follows. 

From this bound, from the pointwise convergence of ip^ io ip m Q and from the fact that 
Fg —)■ F' uniformly on every compact interval included in (—1,1) we infer that 

F^{iPf:) F'(ip) a.e. in Q. (4.106) 


We are now ready to pass to the limit in the variational formulation 04.13p . 04.14p of Problem 
Pe _5 (with (5 > 0 fixed) as e —?■ 0, in order to recover the variational formulation of Problem 
fl4.ip - fl4.7p (still with test functions w e D{A^^‘^) and -0 e H for fl4.ip and fl4.3p . respectively). 
This passage to the limit can be carried out in the same fashion as done in the proof of Lemma 
m (except for the last two terms on the left hand side of fl4.17p which are not present in the 
weak formulation of Problem Pe,^), by employing convergences fl4.96p - fl4.100p . fl4.102p - fl4.105p . 

dnnsi). 

Moreover, in this case, as far as the artihcial term in the momentum balance equation is 
concerned, we observe that its contribution vanishes in the limit. Indeed, fl4.103l) . fl4.104p . 
04.1051) . and fl4.100p entail 


-^p"{ip^)m{ip^){Vip^ ■ V/i,)n„ w^x{t)dt / (^-p'\ip)m{ip)(yip ■ Vfi)u, w^x{t)dt = 0, 


as e —)■ 0, for all test functions w E D{A^/‘^) and y G (7^(0,T), where the last identity is due 
to 04.731) and to the fact that p"(s) = 0 for all s E [—1,1]. Hence, the weak formulation of 
fl4.ip - fl4.7p is obtained in the limit with the function p in place of p, since we have proven that 
\p}\ < 1 and on account of p(s) = p(s) for all s E (—1,1). 

The arguments to prove that u E (^^([O, Tj; that the limit [u,p>] attains the initial 

value [^ 0 ,(^ 0 ] and that the energy inequality 04.161) is satished for almost all t E (0,T) are the 
same as in Step I and we omit the details. 

Finally, we can derive the auxiliary energy inequality fl4.69p satished by the limit [u, <p] 
(with p in place of pf) for all uj E kF^’^(0,T) with oj(T) = 0 and w > 0. Indeed, from Step I we 
can write it hrst for every e—approximate solution [u^, 



-p(po)'“o + + 211 ^ 7 ^ 011 ^) 1 ^( 0 ) 


^p{ip,)ul + E^ip,) + ^11 Vp,f )a;'(r)dr 
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> 


2|| + II Vm((y9e)V/iJ^)a;(r)rfr 


(4.107) 


As is Step I, in order to pass to the limit in the term containing the functional E^{ip^) we need 
to show that 


bj'{T)dT / F^{ip^) / u'{T)dT / F{ip). 


(4.108) 


This convergence can now be established in the following way. First, introduce the function 
dehned by 


G,{s) = F^{s) + 


(4.109) 

with Ooo = ||a||L°°(o) and, observe that owing to fl4.77p . G^ is convex on M. Hence, we have 


Gei^e) < / Ge{^) + / G((v9e)(99e — (^), 


Jfi Jq 

(4.110) 


Introduce the sets := {t G (0,T) : oj'{t) > 0} and := {t G (0,T) : oj'it) < 0}. Then, 
multiply 04.1101) 1 by uj'xi<^ and O4.110p o by where xi% are the characteristic functions of 

the sets /^. Integrating in time between 0 and T and summing the resulting inequalities we 
obtain 

[ u'{T)dT [ Gei^Pe) < [ Uj'{T)dT [ G^{ip) + [ u'{T)dT [ G'((^J (yj, - (^) 


+ / uj'{T)dT / G'(v?)(</?£-<yc) < / uj'{T)dT / G^{ip) 

Ji‘^ Jn Jo Jn 

+ G(^(||G(((^e)||L2(0,r;U) + ||G'((^)||L2(0,r;U))||7’e - 7’||L2(0,r;U), 


(4.111) 


where, in the last inequality we have exploited the bound 04.93p . which in particular implies 
that |G((s)| < |G'(s)| for all s G (—1,1) and for every e G (0,eo]. Therefore, on account of the 
estimate ||F)((<y9e)||2,2(o,T;U) < G, which follows from 04.92p and from the bound ||77elU2(o,T) < G, 
and which implies that ||G((<yCg)||i2(o,T;U) < G and thanks to the strong convergence ^ ip 
in L‘^{Q), we have that the second term on the right hand side of 04.1111) converges to zero as 
e —0 and therefore we deduce 

limsup [ oj'{T)dT [ G^{(p^) < lim f u'{T)dT [ G^{(p) = [ u'{T)dT [ G{(p), (4.112) 


e-i>0 Jo 


e^O 


where the last equality follows from Lebesgue’s theorem (use 04.791) and the fact that |v9| < 1 
a.e. in Q and that F^{s) F{s) pointwise for all s G (—1,1)). On the other hand, thanks to 
Fatou’s lemma and to the pointwise convergence F^{(p^) F{(p), a.e. in Q we also have 


uj'{T)dT / G{(p) < lim inf / uj'{T)dT / G^{(pf). 

n /n /o 


(4.113) 
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From 04.112j) and (14.113p . on account of the definition fl4.109p of we get M.lOSp . 

We can now pass to the limit in fl4.107p in the same way as done in Step I, taking into 
account the strong/weak convergences for and /ie . In particular, the strong convergence 

p{ip^)u‘^ —)■ p{ip)u^ in still holds, and we take into account fl4.104p as well. Therefore, 

we have proven that 04.691) (with p in place of p) is satished and the proof of Lemma O is 
concluded. 

□ 


4.3 Step III. Limit as d —> 0 and end of the proof. 

We now want to pass to the limit in problem 04.ip - 04.7p in order to prove that the original 


problem, i.e. 

{pu)t + div(p'U ® u) — 2div(i/((p)Zi)it) + Vvr + div(ri ® J) = + fi, (4.114) 

div(n) = 0, (4.115) 

+ u ■ V<p = div(m((p)V/i), (4.116) 

p = aip — J*(p + F'{lp), (4.117) 

J := —(3m{Lp)Vp, (4.118) 

dii 

w = 0, 7 — = 0, on dVt, (4.119) 

on 

w(0) = uo, <p(0) = <po, (4.120) 


admits a weak solution. To this aim, we approximate problem 04.114p - 04.12Up by a sequence of 
problems given by fl4.ip - fl4.7p with initial data uq and ipos, where ipos are chosen according 
with the following 

Lemma 6. Given <po £ T°°(fl) with F{ipo) G and |^ol < 1? there exists a sequence 

{(pots) C D{B) with F{(pos) £ T^(fl) and |^ 0 (sl < 1 

<5||V<po<5ir 0 7^05 ^ <Po in H, as 6 ^ 0. (4.121) 

Proof. Take <po 5 ^ F{B) given by 

ipos ■= {I + y/SB)~^ipo. 

Then we have pos + VdBpos = Po and so p^^ + y/dp^s = Po, whence 

l7ol < 1- ('*■122) 

1 + V 0 

Moreover, by introducing the convex function G dehned as in 04.1091) 

G(5) = F{s) + 
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and by multiplying the relation = —VSBifos by G'{(pos) in we obtain 

[ G'{^pos){'^os — V^o) = — f G"f G'{ipos)'^os < 0 , 

Jfi Jq Jq 

owing to the fact that G' is monotone nondecreasing and that it is not restrictive to assume 
that -F'(O) = 0 (and hence G'(0) = 0). Therefore, the convexity of G yields 


G{^pos) < 



G'{^o5)i^os — V^o) < 



(4.123) 


and, since F{(po) G from this last inequality we deduce that F{(pQs) G L^(r 2 ) and hence 

that (fos G with |(^o 5 | < 1 n- 6 . in hi. In order to deduce the hrst of fl4.12ip . notice that 

we have 

Il¥^05|p + ^\Wq6\\1 = (‘d0,‘^05) < 

whence 


<^||V 9 ? 0 (s|P ^ 0, as (5 —)■ 0. 

Finally, the convergence in the second of 04.1211) follows immediately from the theory of maximal 
monotone operators. □ 


Now, choosing the initial data ipos as given by Lemma [HI thanks to Lemma Owe know that 
for every 5 > 0 Problem admits a weak solution [us, ’-ps] having the regularity properties 
04.70p - 04.73p and satisfying the energy inequality 

j ]^p{ips)u] +E{ips)+ ^-\\S7ip5f+ 2 ^ \\^Jv{lps)DusfdT + 5 f dr 

+ WVM^VfisW'^dT < j^^p{(pos)ul +E{(pos)+ ^\\Vipos\\'^ + {h,us)v^ijT, (4.124) 

for almost all t G (0, T). Since <5 is no longer hxed here, the uniform with respect to 6 estimates 
that we shall be able to deduce for the sequence of [u^, will be weaker than the estimates 
obtained in Step I and Step 11. Nevertheless, these estimates will turn out to be enough to pass 
to the limit in P 5 . More precisely, noting that the right hand side of 04.1241) is bounded due 
to Lemma [ 6 ] (see the hrst of 04.1211) and 04.1231) 1. the only uniform in 5 estimates that we can 
write from 04.1241) are now 


p5||L°°(0,T;Gdi„)nL2(0,T;ydi„) < C, 

(4.125) 

p5||i°°(0,T;LP(O))nL2(0,T;y) < C, 

(4.126) 

W'^ P'S\\l^{0,T-,H) < G, 

(4.127) 

IP(7’<5)||l°°(0,T;L1(O)) < G. 

(4.128) 
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Here, we have used (Al)-(A3), the bound < 1 which yields p{(ps) > p* '■= niin(pi,p 2 ) > 0, 
and the fact that, as a consequence of 04.761) (taking e ^ 0) we have 

F{s)>C,\sY-Dj,, Vs e (-1,1), 

with p > 3 hxed arbitrary. 

We now need to control the sequence of averages {ps}- To this aim we hrst consider equation 
04. 3 p written in the form 

V(5 T ^ 

and test it by J\fipg[F'{ips) — F'{ps)) (recall that F'{ps) £ T^(0,T;1/)). Arguing as above we 
still get 

iif'(4,j) • vv^i)ii + viWfsWv) 

^ + 11^5 • VpsWv + 1 ), 

where we have used the bound 5^^‘^\\‘P5\\l°°{o,t-,v) F C, which comes from the energy inequality 
04.1241) . Moreover, on account of the estimate 

\Ws\\l'^{o,t-,v') < C, (4.129) 

which can be obtained by arguing exactly as in Step II by comparison in the weak formulation 
of 04.3p . we are still led to an estimate of the form 

||T'(V5) - F'{‘^i)\\L2(0,T-,H) ^ (4.130) 

With this last estimate available we can now apply the argument devised by Kenmochi et ah 
[IS] to deduce a bound of F'{ps) in T^(0, T; L^(r2)). Let us recall some details of this argument. 
For convenience, here we referee to, e.g., m Proof of Theorem 1]. Introduce first the function 

H{s):=F{s) + ^{s-Sof, Vs G (-1,1), 

where Sq G (—1,1) is such that F'{so) = 0 (cf. (AS)). Hence, owing to (A7) FI' is monotone 
and H'{so) = 0. Then, using the fact that |vo< 5 l < 1 nnd exploiting the argument of Kenmochi 
et ah, the following estimate can be established 

m\\F'{^s)\\Lpn) < [ {^5 (4.131) 
Jn 

where 


ps 1711 , 1712 -Pos}, K{pqs) = {ps +max (|F(| + IF2I) + Ooccr), 

[mi,m2] 

with Ks := max{^Q^ — mi,m2 — niax{so — mi,m2 — Sq}, and mi,m2 G (—1,1) are 

two constants that are independent of <5 and fixed such that mi < Sq < m2 and mi < < m2 
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for all 6 > 0. Due to fl4.122p we see that mi, m 2 can be fixed in this way. Indeed, it is enough 
to £x — 1 < mi < min{ —I^qI, sq} and max{|^Q|, sq} < m 2 < 1. Moreover we can see also that 
the constants ks and hence K(Jpq^) are all uniformly bounded (from below and above) with 
respect to S. In particular, we have ^5 < 1 and min{ —|^ol ~ mi, m 2 — |^ol} < h <5 < 1) for all 6. 
Hence, the constants are bounded by a constant which depends only on Tp^ (and on F, 

J and D). 

By combining fl4.130p (which holds also with H in place of F) with fl4.13ip we deduce the 
desired bound 


and this provides the control ||/i5||L2(o,T) < L, as well as the control 

||h5|U2(0,T;V) < C, (4.132) 

with a constant C now depending also on Tp^, which is derived by using (14.1271) and Poincare- 
Wirtinger inequality. 

Let us now deduce an estimate for the sequence of time derivatives {p{ips)us)t- Recalling 
the weak formulation of (14.11) . for every w E we have 

{{p{ips)us)t,w)^^^,,/^^ = {p{ips)us(^us,Dw) -2{u{ips)Dus,Dw) - 6{A^^‘^us, 

+ / Us - {Js- V)w - ((psVps, w) + (h, (4.133) 

Jn 

where 


Js = -l3m{(ps)Vps- 

From (14.1331) . on account of (14.1251) . (14.1271) and of the bound (5^'^^||'tA5||L2(o. ,t-d(a^/‘^)) < C which 
comes from the energy inequality (14.1241) . it is easy to deduce 

II (p(7’<5)'“(5)t||L2(0,T;D(A3/2)') < C*- (4.134) 

The uniform in 5 bounds (I4.125p . (14.1261) and (14.1321) imply the existence of u, p and p such 
that up to a subsequence we have 


Us ^ 

u. 

weakly* in L°°( 0 , T; Gdiv), 

(4.135) 

Us 

u. 

weakly in L^(0, T; R^j^), 

(4.136) 



weakly* in L°°(0, T; L^’(D)), 

(4.137) 



weakly in L^(0,T; V), 

(4.138) 

Ps 


weakly in L^(0,T; V). 

(4.139) 


Now, we have 

V{p{ps)us) = p{ps)^us + p'{ps)^^s ■ Us, 
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and since us is bounded in L°°(0, T; fl L^(0, T; and hence also in 

then we see that p'{ips)Vcps ■ us is bounded in which implies the bound 

||p((p5)'U5||^5/4(0^T;Wl’5/4(O)3) — (4.140) 

By exploiting the compact embedding and the fact that we have 

2^15/7-7(-^^3 D{A^/‘^y, since D{A^/'^) (here 7 , 7 ' > 0), by Aubin-Lions lemma 

from fl4.134p and fl4.140p we infer that 

p^ip^jUi —)■ p{(p)u, strongly in T; 7 > 0. (4.141) 

Now, observe that 04.1291) and 04.1261) imply that up to a subsequence we have 

ips If, strongly in L‘^{Q), and pointwise a.e. in Q, (4.142) 

p{(ps) —t p{^), strongly in L'^{Q), 2 < Vg < 00 . (4.143) 

Hence, since us ^ u weakly in L^°/^(Q)^, then we have p{p^)us p{'^)u weakly in {QY 

(7 > 0) and by comparison with 04.14ip we get p{(p)u = p{(p)u. Hence, from 04.1341) and 04.1291) 

we have 

{p{ifs)us)t ^ {p{(f)u)t, weakly in L2(0,T;A)(A^/^)'), (4.144) 

if'g ipt, weakly in L^(0,T; V'). (4.145) 

We are now ready to pass to the limit in the variational formulation of Problem P 5 and therefore 
to prove that the original problem 04.114p - 04.120p admits a weak solution. To this goal, observe 
that, due to the interpolation embedding 

L'-(0,T;G,,,) nL2(0,T;Hd*.) U{Q,T -2 < r < cx), 

we have the bound for the sequence of us in L^(0, T; L^‘^/^^(r2)^) and hence 

Us u, weakly in L®(0, T; (4.146) 

The strong and weak convergences 04.1411) and 04.146p . respectively, are enough to pass to the 
limit in the term Ai'v{p{(ps)us ® us) (cf. 04.133p : recall that the test function w E D{A^^‘^) ^ 
and hence Dw E L°°(fl)^^^). 

Furthermore, by employing the pointwise convergence us ^ u a.e. in Q, which follows from 
the strong convergence 04.1411) and from the pointwise convergence 04.142p for ps, we deduce 
the strong convergence 

Us ^ u, strongly in L\0,T-,Gdiv) (4.147) 

(actually the convergence of us to u is strong also in On the other hand, since 

m{ps) —t m{p) strongly in L^{Q) for all g < 00 , by using 04.1391) and the bound || J 5 || 2 , 2 (q )2 < O, 
we obtain 

Js ^ J ■= —(3m{p)Vp, weakly in L‘^{Q)^. (4.148) 
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The strong and weak convergences fl4.147jl and fl4.148jl are now enough to pass to the limit in 
the (variational formulation of the) term div('U 5 0 Jg). 

As far as the term 6A^us is concerned, we observe that from the energy inequality 04.1241) 
we have the bound 5^^'^\\A^^‘^us\\L2(^o,T-,L^{n)^) ^ C and therefore the contribution of this term in 
the variational formulation of vanishes as 5 —)■ 0. Indeed, recalling the standard argument 
to pass to the limit in 04.133p by hrst multiplying it by a test function x G C“(0,T) and then 
integrate the resulting identity in time between 0 and T, we have, for all w G 

fT 

/ {A^/^U5,A^/‘^w)x{t)dt < \/5C'||A3/2^|| ^ 0, as h —)■ 0. 

0 

Regarding the term —6Aips, we see that also the contribution of this term vanishes as 5 —)■ 0. 
Indeed, from 04.124p we have the bound < C*, and therefore for every test 

functions C, and x G C'“(0,T) we have 


(V(^,, VC)x{t)dt < veil ^0, as 5 ^ 0. 


(4.149) 


Now, from the variational formulation of fl4.4p we have 


/ {fJ‘S,0x{t)dt= {aips-J*ips + F'{ips),Ox{t)dt + 6 {V(ps,'^C)x{t)dt, (4.150) 

Jo Jo Jo 

for all e £ R and all x G C^(0)^)- In order to pass to the limit in the term containing F'{ips) 

of 04.1501) we need to show that |(^| < 1 a.e. in Q (observe that up to now, from the pointwise 

convergence 04.1421) and the strict bound \ips\ < 1 for all 5 we only know that |</?| < 1 a.e. in 

Q). To this purpose, we can employ exactly the same argument as in [271 Proof of Theorem 1], 

using the pointwise convergence tps yj, the bound ||RX<^ 5 )||l 1 (q) < C{}Pq), and assumption 

(A8). 

Therefore, since F' is continuous on (—1,1), we have F'{(ps) —t F'{(p) pointwise in Q. This 
pointwise convergence, together with the bound \\F'{(Ps)\\l 2 (^qt.^ < C (which follows from 
fl4.130p and from the bound of 7 I 5 in I^^(0, T)) yield, up to a subsequence 


F'{ips) F'{ip), weakly in L^{0,T-, H). 


(4.151) 


We can then pass to the limit as 5 —)■ 0 in fl4.150p . using fl4.149p and the weak convergences 
fl4.15ip . (14.1391) This shows that fl4.117p is satished by the limit functions and /r. 

The argument for passing to the limit in the other terms of Problem P 5 is straightforward 
and therefore, by letting 5 —)■ 0 from the weak formulation of Problem P^ we hnally recover the 
weak formulation of the original problem fl4.114l) - fl4.120p with test functions w G D{A^F'j 
■0 e R for fl4.114p and fl4.116p . respectively. 

We can now observe that, by density, the weak formulation of (14.1141) holds also for every 
test function w G D{A). Moreover, for all w G D{A) we have 


(u I 


J) : Dw 


< C'll^llL3(o)||J||Nlb(A) < C|| 


u 


11/21 


u 


11/2 I 
lL6(n)3| 


\w 


llD(n), 
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whence we obtain 




div(n (g) J) e T; D{A)'). 


Indeed, we have 


(pu 0 u) : Dw 




for all w G which entails 

div(p'u 0 w) e L^/^(0,T; I/jiJ. 

It is also immediate to check that the terms —2div(z/(9?)Zi)n) and pVip in 04.114p are in 
L^(0,T; Hence, fl3.6p is satished for all w G D{A) and we have 

{pu\EL^/HO,T;D{Ay). 


Let us now prove the weak continuity u G Cw{[0,T];Gciiv)- To this purpose note that from 
(14.125p and fl4.134p we have the boundedness of 

p{(ps)us in L°°(0,T; L^(r2)^) and 

p{ips)us in H\0,T-D{A^/^y) -G C{[0,T]-D{A^/^y). 

Therefore, Lemma [2] and fl4.14ip imply that p{(p)u G C^([0, T]; L^(r2)^) and this leads to the 
weak continuity of u in Gdiv as in Step I, on account of the strong continuity ip G C'([0, T]; H). 

The argument to prove that u and p attain the initial values Uq and poj respectively, follows 
the same lines as in Step I and we omit it. 

Finally we prove the energy inequality fl3.9p . First, we know from Step II that [u^^ps] 
satishes fl4.69p . namely 

( [ \pi^os)ul + E{pos) + ^||V</?o5ir)a;(0) + f f ^p{p5)uj + Eips) + ^\\Vps\y^uj'iT)dT 

T T 

> j (^2||a/^<^L>W 5|P + + ||y/m(^V/i5||^)a;(r)dr - y {h,us)vai^uj{T)dT, 

(4.152) 

for every co G kF^’^(0,T), with u){T) = 0 and ca > 0. Again, to pass to the limit in 04.152p we 
need to show that 

u'{T)dT I E{p). (4.153) 

Jn 
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This convergence can be established by the same argument as at the end of Step II, by exploiting 
the fact that F is a quadratic perturbation of a convex function and by using the bound of 
F\ip^) in L^(0,T;Tr). Note that if F is bounded, then fl4.153p follows at once by directly 
applying Lebesgue’s theorem. Concerning the other terms in fl4.152l) . we have 

[ u'{T)dT [ ^p{Lps)uj [ u'{T)dT [ (4.154) 

Jo Jn Jo Jn ^ 

thanks to the strong convergences —)■ n in L^{QY (which follows from the bound of 

in and pointwise convergence) and p{i^s) p(</?) in L^{Q) for all g < cxo (take 

q = 3). Convergence fl4.154p can also be justihed by observing that we have also ^p((ps)u^ —t 
In almost all r G (0,T) and by applying Lebesgue’s theorem. 

Furthermore, we have 


/ / -p((po)ul 

thanks to Lebesgue’s theorem, 

[ F{^pos) < f F((po) + ^(11(^011^ — ll^^o^lP) ^ [ F{ipo), as (5 —)■ 0 , 

Jn Jn ^ Jn 

thanks to fl4.123p and to the second of fl4.12ip . and 

I < C'a;||<^5|lL2(0,T;V) < aS (5 0, 

0 2! 

thanks to the bound fl4.126p . By also using the hrst of 04.1211) and weak lower semicontinuity 
of norms (the second term on the right hand side of the inequality is simply neglected) we can 
now pass to the limit in 04.1521) and obtain the following integral inequality 


£^(0)ci;(0) + / 8{t)ijj'{T) dT > / V{T)dT, 


(4.155) 


satished for every u G hF^’^(0,T), with u)(T) = 0 and a; > 0, where the functions S and V are 
given by 

I ^p{ip)u^ + E{ip), V{t) := 2\\^/^J{^Duf + \ym{p:)V pf - 


for all t G [0,T] and almost all t G (0,T), respectively (for simplicity of notation, we omit the 
indication of time t on the right hand side). Let us check that S = (p(-)) : [0, oo) —)■ M is 

lower semicontinuous. Indeed, we know that p{(p)u G Ci„([0, T]; L^(r2)^) and moreover we have 
p{ip)~^F £ C([0,T];Fr). Therefore p(ip)u = p((p)“^/^p((p)it G Cu,([0, T]; L^(r2)^), on account 
of the boundedness of p and this proves the lower semicontinuity in time of the hrst term of 
£. The lower semicontinuity of E{ip{-)) : [0, cxo) —)■ M is a consequence of the fact that F is a 
quadratic perturbation of a convex function (see [2Sl Lemma 2]). 
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The energy inequality fl3.9j) now follows by applying Lemma [3] to fl4.155jl . 

The proof of Theorem [T] is now complete. 

□ 

Remark 5. The two dimensional case. The proof of Theorem [T] can be obviously carried 
out in two space dimensions as well. In particular, the strong convergences that can be obtained 
for the sequences of approximate solutions at each step by using interpolation and Aubin-Lions 
lemma will generally hold in stronger norms with respect to the 3D case. Nevertheless, after 
passing to the limit as 5 —?■ 0 the weak solution we get for system fll.ll) - fll.8p still has no more 
than the regularity given by fl3.1l) - fl3.3p . fl3.5p and by the second of fl3.4p . Only the regularity for 
{pu)t improves a bit. Indeed, by a comparison in the weak formulation fl3.6p with test function 
w G for 1 < s <2, using interpolation and Gagliardo-Nirenberg inequality in 2D it is 

not difficult to see that 

{pu)t e L2 -t'( 0,T;D(A)') nL2/(3-")(0,T;D(A"/2)'), 

for every 0 < 7 < 1 and every 1 < r < 2. Notice that we cannot set r = 1. Indeed, in the weak 
formulation 03.61) we cannot take the test function w in Vdiv due to the extra-term div(n G) J) 
and to the fact that u does not belong to Therefore, even if the regularity for {pu)t 

slightly improves in 2D, this is not enough to show the validity of the energy identity. For the 
same motivations also uniqueness of weak solutions in 2D is not known. 

In conclusion, in order to get an improvement of the results in 2D concerning uniqueness and 
validity of the energy identity, we hrst need to prove a regularity result for system fll.ip - fll.8p . 
namely, to establish existence of solutions that are more regular than the ones constructed in 
Theorem [1] (a regularity assumption on the initial data will then be required). This will be the 
subject of a future contribution. 
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